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Abstract. The impact of current on static and kinetic depinning fields of a domain wall in an one 
dimensional ferromagnetic nanostrip is investigated by solving the Landau-Lifshitz-Gilbert equa¬ 
tion with adiabatic and non-adiabatic spin-transfer torques analytically and numerically. The results 
show that in the absence of current, the static depinning field is greater than the kinetic depinning 
field and both the depinning fields decrease by the increase of current applied in a direction opposite 
to the direction of the applied field. Both the depinning fields can also be tuned by the current to 
make them equal. 
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1. Introduction 


In the recent years, research in the motion of domain wall in ferromagnetic nanostrip has 
been focused more due to its practical applications such as magnetic logic III] and memory 
devices|2]. The manipulation of domain wall in r nag netic nanostrip can be achieved by 
the application of magnetic field|3] and/or currentUl^j. By applying the magnetic field 
and current, the domain wall moves rigidly up to the critical value known as Walker limit. 
Above the Walker limit, the motion of the domain wall is not regular and exhibits oscil¬ 
latory behaviorjl h l Ml . Conventionally, the position of domain wall can be controlled by 
ion irradiation U[lll], nearby nano-magnets ll3 - [l5ll and introducing artificial geometrical 
constraints such as notches in the long edge of the wire] 3-^. The displacement of a 
domain wall from one notch to another notch can be achieved by depinning the trapped 
wall from one notch and pinning the moving wall in the other notch and its pinning mech¬ 
anism could be understood by static and kinetic pinning respectively. The trapped domain 
wall can be depinned from a notch when the applied magnetic field is above the threshold 
magnetic field, which is called the static depinning field. Similarly, the moving domain 
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wall cannot be trapped, when it crosses the notch if the applied magnetic field is above the 
threshold magnetic field, which is known as the kinetic depinning field. The static depin- 
ning field and kinetic depinning field are independent of the shape of the notch whether 
it is rectangle or triable in shape and the kinetic depinning field depends on the chirality 
of the domain walll 17 .B. 

Instead of producing a notch by single triangle, one can create the notch by two sym¬ 
metrical triangles in the nanostrip which has more advantage than the single one because 
of the irrelevance of static depinning field on the chirality and propagation direction of a 
domain wall lfB . Recently, Sung-Min Ahn et 0^12^ 23] have found the following results 
for the notch created by two symmetrical triangles. The static and kinetic depinning fields 
decrease with the increase of the width of the nanowire ll22ll and the static depinning field 
is greater than the kinetic depinning field JB . For the above results, they have modelled 
the pinning field created by the notch as a step function. Though there are plenty of stud¬ 
ies in field driven domain wall pinning, the study of domain wall pinning in the presence 
of current is limited. When the applied current is above the threshold, the trapped domain 
wall in a notch can be untrapped, and the corresponding current is known as static depin¬ 
ning current. It increases when the notch depth ll^ and notch anglejB are increased. 
And also, the static depinning field reduces when the current is applied in a direction 
opposite to the direction of the domain wall propagation ll^ . The static depinning cur¬ 
rent increases or decreases with the increase of doping concentration or magnetic field 
respectively 

In this paper, the dynamical equation of domain wall with a symmetrical notch has been 
solved analytically and numerically. The results show the variation of static and kinetic 
depinning fields and the equivalence of both depinning fields in the presence of current. 
The paper is organised as follows: The model of the domain wall along with the pinning 
field is discussed in Section 2 and the corresponding equations for the velocity, width and 
excitation angle of the domain wall are derived analytically in Section 3 in the presence 
of current along with the pinning field. The numerical results are explained in Section 4 
and finally the results are summarized in Section 5. 


2. Model 



Figure 1. (a) A Schematic diagram of the ferromagnetic nanostrip having an artificial 
notch in the form of two symmetrical triangles which is taken as our model and its mag¬ 
netization varies only along x-direction. The current density J and external magnetic 
field H have been applied in negative x-direction and positive x-direction respectively, 
(b) Representation of magnetization in terms of spherical coordinates. 

We consider an infinitly long ferromagnetic nanostrip with a single Neel-type trans¬ 
verse domain wall as shown in FIG. 1(a) as our model to study the phenomena of static 
and kinetic pinning of the domain wall. The pinning field is introduced by etching the cen¬ 
ter of the nanostrip with an artificial notch which is in the form of symmetrical triangles 
as shown in FIG. 1(a). This symmetrical nature of the notch is more advantageous than the 
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asymmetrical notch because of its ability to create an unique static depinning field irre¬ 
spective of the domain wall polarity and its propagation direction. The two domains of the 
ferromagnetic nanostrip are separated by a domain wall with opposite magnetizations(M) 
along the easy axis(x-axis) and the arrows in the strip show the direction of magnetiza¬ 
tion. The external time varying magnetic field(H(t)) and current(J(f)) are applied along 
the positive and negative x-direction respectively. In FIG. 1(b), x, y, z represent the unit 
vectors along x,y and z directions respectively and the angles 6 and $ refer to the angle 
between the magnetization vector M and the positive x-direction and the angle between 
the projection of M in yz-plane and positive y-direction respectively. Assume that the 
variation of M takes place only along the x-direction and the angles 0 and T* imply the 
deviation of the magnetization vector from positive x-direction and xy-plane respectively. 
The dynamics of domain wall is understood through the spatial and time variation of 
magnetization governed by the Landau-Lifshitz-Gilbert (LLG) equation in the presence 
of field and cuiTent. The LLG equation is written asfsll 


i9M(a;, t) 

m 


= -7M X He// ■ 


a 5M 

—M X 

Me dt 


Ht) 

M? 


M X M X 


5M 

dx 


5M 

—M X , 
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M = (M„ My, M,)- |M|2 = Ml + Ml + Ml = Ml, 


(la) 

(lb) 


where, 


h{t) = 


eMs 


c{t) = 


(2a) 

(2b) 


Here, 7 is the gyromagnetic ratio, a is the Gilbert damping parameter, Ms is the satu¬ 
rated magnetization. He// is the effective field. b{t) and c{t) represent the magnitude of 
adiabatic and non-adiabatic spin-transfer torques respectively, which include the interac¬ 
tion of conduction electrons and the local magnetization. P is the polarization, J is the 
magnitude of current density, hb is the Bohr magneton, e is the charge of electron and 
a- 0.01)11 is the ratio between c and b. The adiabatic spin-transfer torque is corre¬ 
sponding to the reaction torque on the magnetization produced by the spatial variation of 
the spin current density And the non-adiabatic spin-transfer torque corresponds to 
the reaction torque on the magnetization due to the continuous space variation of spatially 
mistraking spins between conduction electrons and local magnetizationH. The initial and 
final velocity of a domain wall is controlled by adiabatic and non-adiabatic spin-transfer 
torques respectively fll. The effective field He// includes the fields due to exchange 
energy, easy axis anisotropy, external field, demagnetization field and the pinning field 
produced by the notch lEl is given by 

o A / fT \ 

= + 

where. 


U{x) = 1, for 0 < X < 5, (4a) 

= 0, otherwise. (4b) 

Here, A represents the exchange interaction coefficient, Hk represents the magnetocrys¬ 
talline anisotropy coefficient, H is magnitude of external magnetic field and AttM^ is the 
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demagnetization field. Hp is the pinning field, U(x) is the step function used to introduce 
the pinning held in the region 0 < x < S and S is the range of the pinning held. 


3. Analytical solutions for domain wall parameters 

In order to understand the dynamics of the domain wall in pinning held, we try to solve 
Eg. (flat. As Eg. ([Tall is a highly nontrivial vector nonlinear evolution equation, it may be 
difficult to solve the same in its present form. Hence, we rewrite Eq. ( [Tat in terms of the 
polar coordinates by using the transformations = Ms cos 0, My = Ms sin 0 cos <i>, M^ = 
Ms sin 9 sin $ and the resultant equations read 


89 . 2jA 

— +asm9-—= —— 
dt dt Ms 


„ ^89 8^ . 

2cos6»—— +sin6»-^ 
OX ox ox^ ) 


do 

—27r7Ms sin0 sin2$ + hit )——h cit) sin0-—, 

8x 8x 


89 . 


_ 27A / 8^9 

Ms i 8x‘^ 


5$' 

— sin 9 cos 9 ( 

8x , 


(5a) 


— 'y[H{t) — HpU{x)] sin0 — ^ + AttMs sin^ $) sin 20 

89 

- 6(f) sin 0—+c(f) — . (5b) 


Eqs.Q can be solved by using the following trial functions introduced by Schryer and 
WalkerHU. 


9{x, t) = 2 tan ^ exp 

^{x,t) = (6b) 


where, X (f) is the position of the center of the domain wall and W (t) is the width of the 
domain wall. Eo.d^ assumes that as time goes on, the domain wall moves without any 
change in the static prohle(the spatial variation of 0 in the absence of current and held) 
and only with the change in its width. And Eo. (l6bl) assumes that the domain wall excites 
from the xy-plane with time and the excitation is independent of space. (j){t) is an angle 
between the projection of magnetization(yz-plane) and positive y-direction, which can be 
called as excitation angle. 

On substituting ^ from Eo. dSbl i in Eo.dSall. we get 


(1 + a )sin0— = 


8^ 2 A 7 
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a Sint 
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2a cos 0 
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9$ do 

+ (1 + a^)b{t) sin0 — + (a - ^)b{t) 
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Similarly, by substituting ^ from Ea.(l5all in Ea. (l5bl i. we get 


2. do 2A'y 


+ “ - M. 


. , 92 $ ^ . 50 9 $ d^e . . , / 9 $ 

sin 0 + 2 cos 0 — ——h a-—^ — asmO cos 0 | 

ax ax ox^ 


dx 


— 2771 Mg sin 2$ sin 9 — ^ [Hk + AttMs sin^ $] sin 29 — a'y[H(t) — HpU (a;)] sin 9 


do 

- (a - ^)b{t) sin 0 — + (1 + af)b{t) 

Erom Eas.d^ and dbbl) . one can derive the following identities at x = X. 
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(9a) 

(9b) 

(9c) 

(9d) 

(9e) 

(9f) 

(9g) 

(9h) 


The reduced form of Eq.O is obtained by substituting x = X(t) in it. Eurther, it can be 
simplified by using the identities given in Eqs.® and we can obtain, 


(1 + - HpU{X) - 2naM, sin 20) + ■ 

( 10 ) 


Similarly from Eq.® we can derive 

= ^ = 77 ^ 7 ^%( 27 rMg sin 20 + a[iT(f) - HpU{X)\) - b{t) 
dt (1 + a^) 


^ + ^2^ ■ 

(11) 


The width of the domain wall is obtained by differentiating Eq. O with respect to x and 
reducing it at x = X (t) and the reduced equation can be simplified using the identities 
given in Eqs.®. 

Wit) = VEo [1 + (4^Mg/Mfe) sin^ 0] . (12) 

Here, Wq = y/2A/HkMs is the initial width of the domain wall. Eas. lfTOl) . (fTTI) and 
(fT2l l give the variation of the excitation angle, velocity and width with respect to time re¬ 
spectively in the presence of current and external field as time passes. Since Eqs. (fTTTl and 
dEll involve 0, in order to find the velocity and width of the domain wall, it is needed to 
solve Eg. ([Tol l. As Ea. tfTOl i is a highly nontrivial nonlinear evolution equation it is difficult 
to solve the same analytically. Hence, we solve Ea. ifTOl i numerically and the results are 
discussed in the forthcoming sections. 
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4. Numerical Results 

The pinning and depinning of the domain wall can be understood numerically by inte¬ 
grating the dynamical equation (fTOl l for the excitation angle (p by using Runge-Kutta-4 
algorithm with the initial conditions cj){0) = 0 and using the experimentally measured 
values of the material parameters for Cobalt nanostripes as given by Mg = 14.46 x 
Am-^, Ms = 1.8 X 10^/471 Oe, A = 2x10-^^ Jm-\ j = 1.9x10^ Oe-^s-\ Hk 
500 Oe and P = O.35|01- The time varying field and current density are taken in the fol¬ 
lowing form. 


H{t) = [H'/h)t for 0 < t < h, 
= H' when t > h, 
J{t) = {J'/h)t for 0 < t < h, 
= J' when t > h, 
h = 10-12. 


(13a) 

(13b) 

(13c) 

(13d) 


Eas. lfOl) represent the external field H and the current density J in the form of linearly 
increasing pulse with the duration of 1 ps and after that both are maintained as constant. 
H' and J' can be referred as the saturated external field and current density respectively. 
If the value of H' or J' is just above the static depinning value, the wall is depinned 
from the pinning field region and travels along the direction of the field (or) along the 
direction opposite to the direction of the current. To avoid the non-zero initial velocity of 
the domain wall when it enters the pinning field region to find the static depinning field, 
the values of H and J are increased from zero. Otherwise, there will be no difference 
between the kinetic and static depinning. h decides the rate of change of the field and 
current which has been fixed as 1 ps throughout this paper. In order to understand the 
static and the kinetic depinning of field and current, the displacement X (t) of the domain 
wall is obtained numerically by applying Simpsons’s 3/8 rule as follows. 



(14) 


4.1 Static pinning and depinning under field and current: 

In this section, the details of the systematic investigation of the static pinning and de¬ 
pinning of a domain wall in the presence of either field or current are discussed. Ini¬ 
tially, the domain wall is placed at a position a; = 0, just before the pinning field region 
(0 < a; < i5) with a pinning field of strength Hp=lQQ Oe. The domain wall is moved from 
rest by applying a time varying field 7T(saturated field; H') along the positive x-direction 
or by applying a current density J(saturated current density: J') along the negative x- 
direction. To understand whether the wall moves beyond the pinning field region or not 
under field(current), the displacement of the domain wall is numerically calculated from 
Ea.(fT4t for If the strength of the saturated field H' (or) current density J' is 

not sufficient to move the domain wall beyond the pinning field region, then the wall is 
trapped or pinned within the pinning field region created by the notch in the nanostrip. 

The displacement of the wall X{t) is obtained numerically from Ea. (fT4b for different 
strengths of saturated external field namely H'= 99.9 Oe, 100.0 Oe and 100.1 Oe for a 
pinning field strength T/p= 100.0 Oe and for a pinning field range of S=Wo=23.52 nm is 
shown in EIG.2(a). When H' <= 100.0 Oe, the domain wall is pinned in the pinning 
field region and the corresponding displacement is constant. Once the field H' is just 
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above the 100.0 Oe, the wall crosses the pinning field region and displaces linearly as 
time passes on as shown in FIG.2(a), which indicates that the static depinning field(iFsdp) 
of the wall driven by field in the absence of current is 100.0 Oe and it is equal to the 
given pinning field strength. In a similar way, the static depinning current density) Jg^p) 
is understood from FIG.2(b), where the plots showing the displacement of the domain 
wall against time have been plotted for different values of saturated current density: J' = 
-58.7x 10® A/cm^, -58.8 x 10® A/cm^ and -58.9x 10® A/cm^ for Hp=\00 Oe and 5= 23.52 
nm. For the saturated current densities -58.7 x 10® A/cm^ and -58.8 x 10® A/cm^, the wall 
is pinned in the pinning field region, whereas for J'=-58.9x 10® A/cm^ the wall moves 
out of the pinning field region as shown in FIG.2(b). Hence, the static depinning current 
density Jgdp of the wall driven by current when the field is absent is given by -58.8 x 10® 
A/cm^. 



Figure 2. (color online). Displacement X{f) of the domain wall against time for (a) 
field driven and (b) current driven cases. The pinning field strength is fixed as Hp=\QQ 
Oe and the range of the pinning field as 5 = lTo=23.52 nm. 


4.2 Kinetic pinning and depinning under field and current: 

In the case of static pinning and depinning, the domain wall is moved from rest into the 
pinning field region whereas the kinetic pinning and depinning of the field and current are 
studied by allowing the moving domain wall into the pinning field region. Therefore, the 
initial position of the wall is taken far away from the pinning field region along negative 
x-direction in such a way that the domain wall moves with constant velocity when it enters 
into the pinning field region. The strength of pinning field and its range is kept as same 
as the static case. By applying the field(current) with their saturated values of H'{ J'), the 
domain wall starts to move and reaches the pinning field region with constant velocity. 
The wall would not be pinned if H'(J') is just above a threshold value, otherwise it would 
be pinned. This threshold value of field and current are referred to kinetic depinning 
fie\d(Hkdp) and kinetic depinning current density) J^^p) respectively. 

The displacement(Ai) of the domain wall versus time has been plotted for the field 
driven and current driven cases, are shown in FIGs.3(a) and 3(b) respectively. The initial 
position of the domain wall is set as Ar(0)=-1000 nm and the strength and range of the 
pinning field is considered as 100 Oe and S — 23.52 nm respectively. In the absence 
of current, when the saturated field is equal to or below 36.7 Oe, the wall moves from 
x=-1000 nm and settles at 0 nm where the notch(pinning field region) is located. When 
the saturated field H' is increased to 37.8 Oe, the domain wall moves beyond the pinning 
field region as shown in FIG.3(a). In a similar way, in the absence of the field, when the 
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Figure 3. (color online). Displacement of the wall against time for (a) field and (b) 
current driven cases. The velocity of the wall with respect to the displacement (c) at 
the depinning field, (d) above the depinning field for the field driven case and (e) at the 
depinning current density, and (f) above the depinning current density for the current 
driven case. The strength of the pinning field and the range are set to be Hp =100 Oe 
and 5 = 23.52 nm for all the plots. 


saturated current density is equal to or below the value of -78.9x10® A/cm^, the displace¬ 
ment of the wall X(t) settles in the notch at 0 nm. When the saturated current density 
J' is increased to -79.0x10® A/cm^, the domain wall moves beyond the pinning field re¬ 
gion as shown in FIG.3(b). From both the figures (3a and 3b), one can observe that the 
kinetic depinning field H^dp and kinetic depinning current density J^dp are 36.7 Oe and - 
78.9 X10® A/cm^ respectively. It is interesting to observe from the static and kinetic cases, 
in the absence of current, the static depinning field is greater than the kinetic depinning 
field whereas in the absence of field the static depinning current density is lower than the 
kinetic depinning current density. 

The pinning and depinning of the moving domain wall in the pinning field region can 
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be observed from figures 3(c)-(f) plotted between the velocity and the displacement of 
the domain wall. FIG.3(c) and 3(d) show the plots corresponding to the saturated fields 
H' = Hkdp and H' > Hkdp in the absence of current respectively. When H' is equal to 
the kinetic depinning field(36.7 Oe), the velocity of the wall attains the constant value be¬ 
fore reaching the pinning field region 0 < x < 23.52 nm and then it drastically drops and 
oscillates in the pinning field region and finally, the velocity reaches zero at x=0 nm and 
the domain wall comes to rest as shown in FIG.3(c). When the field H' is increased from 
36.7 Oe to 36.8 Oe, the velocity of the wall drops quickly when it enters into the pinning 
field region and regains the same velocity when it comes out of the pinning field region, 
which is observed from FIG.3(d). Similarly, the velocity of the domain wall against dis¬ 
placement is plotted for the saturated current densities equal to and above the kinetic 
depinning current density( Jfc£ip=-78.9 xlO® A/cm^) in the absence of field as shown in 
FIGs.3(e) and 3(f) respectively. FIG.3(e) shows that when J' — Jkdp the velocity of the 
domain wall decreases to zero when it enters into the pinning field region 0 < x < 23.52 
nm. When J' is just above the kinetic depinning current density (-79.0 x 10® A/cm^), 
the velocity of the domain wall decreases towards zero but it does not reach zero in the 
pinning field region and it retains back its velocity after crossing the pinning field region 
which is shown in FIG.3(f). The velocity of the domain wall after crossing the pinning 
field region is found as 700.8 m/s for i7'=36.8 Oe which is slightly above the kinetic 
depinning field whereas in the current driven case the velocity is found as 553.0 m/s for 
J'=-79.0 X10® A/cm^ which is slightly above the value of the kinetic depinning current 
density. The velocity of the domain wall to cross the notch is less in the current driven 
case compared to the field driven case. 

4.3 Effect of current on the static and kinetic depinning field: 

The variation of the kinetic depinning field H^dp and the static depinning field Hgdp with 
respect to the saturated current density J' is discussed here. In order to understand the 
static and kinetic depinning fields numerically, the initial position of the domain wall is 
fixed at 0 nm and -1000 nm respectively for a fixed saturated current density. The saturated 
field H' is increased from 0 by an increment of 0.1 Oe. For each and every increment of 
H' Ea. (fT4l) is solved numerically and the values of static and kinetic depinning fields are 
obtained. This can also be repeated by increasing J' from 0 to -80 x 10® A/cm^ with an 
increment of -0.5 x 10® A/cm^ and the results have been plotted in FIG.4. For both the 
cases of static and kinetic pinning, Hp and the range 5 are fixed as 100 Oe and 23.52 
nm respectively. In the absence of current, the kinetic and static depinning fields are 
found to be 36.7 Oe and 100.0 Oe respectively. This proves that the magnitude of the 
static depinning field is greater than the kinetic depinning field which has been evidenced 
in the earlier work lEill . When the saturated current density is increased in the opposite 
direction of the external field, both the depinning fields decrease, especially the static 
depinning field decreases quickly than the kinetic depinning field. At a particular value of 
the saturated current density(=-54.47 x 10® A/cm^), both the static and kinetic depinning 
fields get equal value (16.7 Oe) and above this value of J', the static depinning field takes 
the lower value than the kinetic depinning field. The decrease in static depinning field 
with the increase of | J' | when current is applied in the opposite direction of field implies 
that a pinned domain wall can be depinned with the lower external field in the presence 
of current. Similarly, the decrease in kinetic depinning field implies that the moving wall 
can cross the notch with smaller external field in the presence of current. From FIG.4, one 
can understand that the increase in saturated current density decreases both the depinning 
fields, and also makes the kinetic depinning field to be greater than the static depinning 
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Figure 4. (color online). The static depinning field(red) and the kinetic depinning 
field(blue) with respect to the saturated current density(J') with JTp=100 Oe and 
(5 = 23.52 nm. 


field. For example, when J'=-52.0x 10® A/cm^, the static depinning field is 25.2 Oe and 
the kinetic depinning field is 17.8 Oe, which implies that the static depinning field is 
greater than the kinetic depinning field. However, when J'=-58.0x 10® A/cm^, the static 
and kinetic depinning fields are 3.5 Oe and 15.0 Oe respectively, which implies that the 
static depinning field is smaller than the kinetic depinning field. This shows that, one can 
tune the values of static and kinetic depinning fields and make them equal by tuning the 
current density. Also it can be verified that the static and kinetic depinning currents are 
-58.8 X10® A/cm^(label a) and -78.9x 10® A/cm^(label b) respectively in the absence of 
field as shown in FIG.4. 

5. Conclusions 

In this paper, the pinning and depinning of a Neel-type transverse domain wall in the 
pinning field region or notch created by two symmetrical triangles are reported in the 
presence of field and current. The corresponding LLG equation with the adiabatic and 
non-adiabatic spin-transfer torques along with the effective field is solved analytically to 
derive the velocity, width and excitation angle of the domain wall. The displacement 
of the domain wall is obtained by numerically integrating the velocity equation using 
Sympson’s 3/8 rule. The static and kinetic depinning fields in the absence of current have 
been found as 100 Oe and 36.7 Oe respectively for a pinning field strength 100 Oe in the 
range 0 < a; < 23.52 nm . Also, the static and kinetic depinning current densities have 
been found as -58.8 x 10® A/cm^ and -78.9 x 10® A/cm^ respectively for the same pinning 
field strength and range in the absence of the field. For the field driven case, the static 
depinning field is greater than the kinetic depinning field whereas for the current driven 
case, the static depinning current density is smaller than the kinetic depinning current 
density. The velocity of the domain wall to cross the notch is less in the current driven 
case(553.0 m/s) when compared to the field driven case(700.8 m/s). At a particular value 
of the saturated current density(-54.5 x 10® A/cm^) both the depinning fields are equal and 
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below -54.5 x 10® A/cm^, the static depinning field is greater than kinetic depinning field 
whereas above -54.5x10® A/cm^, the static depinning field is smaller than the kinetic 
depinning field. This work helps to improve the mechanism to control the motion of the 
domain wall by geometrical notches in ferromagnetic nanostrips. 
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